/*******************************************************************************
* Draw map
*******************************************************************************/

use "$edited/final_data", clear
keep if date==22500

xtset, clear

format sh_dose_1_20211214 %12.1f
grmap sh_dose_1_20211214, clnumber(10) fcolor(Blues2) ///
	polygon(data("$edited/state_geo_shp.dta") osize(medthick)) ///
	point(data("$edited/geopoints.dta") xcoord(_CX) ycoord(_CY) select(keep if any_antivax_gps==1) ///
	shape(diamond) fcolor(red) size(*0.3) legenda(on) leglabel("At least one vaccine-skeptic GP")) ///
	legorder(lohi) legstyle(2) legtitle("Population vaccinated with a least one dose" "on 2021-12-14 in % (deciles)") legend(position(11)) name(unconditional, replace)

* pdf version
graph export "$output/map_vaccination_1dose.pdf", replace

* jpg version
graph export "$output/map_vaccination_1dose.jpg", replace	

